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Abstract. If a neutral inclusion is inserted in a matrix contain- 
ing a uniform applied electric field, it does not disturb the field 
r*j | outside the inclusion. The well known Hashin coated sphere is an 

O j example of a neutral coated inclusion. In this paper, we consider 

the problem of constructing neutral inclusions from nonlinear ma- 
terials. In particular, we discuss assemblages of coated spheres and 
the two-dimensional analogous problem of assemblages of coated 
disks. 



1. Introduction 

A neutral inclusion, when inserted in a matrix containing a uniform 
applied electric field, does not disturb the field outside the inclusion. 
The problem of finding neutral inclusions goes back to 1953 when Mans- 
field found that certain reinforced holes, which he called "neutral holes", 
could be cut out of a uniformly stressed plate without disturbing the 
surrounding stress field in the plate |Man53j . The analogous problem 
of a "neutral elastic inhomogeneity" in which the introduction of the 
inhomogeneity into an elastic body (of a different material), does not 
disturb the original stress field in the uncut body, was first studied by 
Ru } Ru98] . 

The well known Hashin coated sphere is an example of a neutral 
coated inclusion [Has6_ 2] . In |HS62a[ IHS62b] . Hashin and Shtrikman 
found an exact expression for the effective conductivity of the coated 
sphere assemblage. For the case of imperfect interfaces, neutral spher- 
ical inclusions have been studied in |LV95] , |LV96a] , |LV96b] , |Lip97a| , 
|Lip97b| , and [IT99J . Examples of neutral inclusions of arbitrary shape 
have been studied in [MSOlj . For other references see |MV10] and 
Chapter 7 of |Mil02j . 



Key words and phrases, neutral inclusions; inclusion assemblages; p-Laplacian. 

1 



2 SILVIA JIMENEZ, BOGDAN VERNESCU, AND WILLIAM SANGUINET 



To begin, we consider a particular coated sphere (see Fig. ([[])), of 
phase 1 in the core and phase 2 in the coating, of core radius r c and 
exterior radius r e , 1 < r c < r e , subject to linear boundary conditions, 
i.e. we apply the linear field E • x = Ex%, (where E = Ee 1 , with 
e 1 = (1,0,0) and x = (xi, rr 2 , ^3)) as a boundary condition to the 
exterior boundary of the sphere, to find a solution u that solves 

aiA p u = (7i V ■ (|Vm| p ~ 2 V«) = (nonlinear) in the core, ^ 
a 2 Au = (T2V ■ Vm = (linear) in the coating 

where the conductivities Oi for i = 1,2 satisfy 00 > o\ > cr 2 > 
and A p represents the p-Laplacian (p > 1), together with the usual 
compatibility conditions at the interfaces (continuity of the electric 
potential and continuity of the normal component of the current). The 
ratio of the inner sphere radius to the outer sphere radius is fixed, 
therefore the volume fractions of both materials 81 and 62 are fixed: 

4/37rr 3 r 3 

0! = -^— f = -|, and 02 = 1-6,. 1.2 




Figure 1. Cross section of a coated sphere 

Our calculations in Section [2] show that we can replace the coated 
sphere with a sphere composed only of linear material of conductivity 
a* (see (12.151) and Fig. [2]). Since the equations for conductivity are local 
equations, one could continue to add similar coated sphere of various 
sizes without disturbing the prescribed uniform applied field surround- 
ing the inclusions. In fact, we can fill the entire space (aside from a 
set of measure zero) with a periodic assemblage of these coated spheres 
by adding coated spheres of various sizes ranging to the infinitesimal 
and it is assumed that they do not overlap the boundary of the unit 
cell of periodicity (see Fig The spheres can be of any size, but 

the volume fraction of the core and the coating layer is the same for 
all spheres. While adding the coated spheres, the flux of current and 
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electrical potential at the boundary of the unit cell remains unaltered. 
Therefore, the effective conductivity does not change. 




FIGURE 2. When a* is chosen appropriately, one can in- 
sert a coated sphere into the medium without disturbing 
the surrounding uniform current field 

This configuration of nonlinear materials dissipates energy the same 
as a linear material with thermal conductivity <t* (see (I2.15P ). 

This paper is structured as follows: Section [2] provides the state- 
ment of the problem and the main result for an assemblage of coated 
spheres, Section [3] extends the results for an assemblage of coated disks, 
Section 0] contains a discussions about the effective conductivity of the 
assemblage and its relation to p, and Section [5] contains the conclusions 
and a description of future work. 

2. Assemblage of Coated Spheres: Statement of the 
Problem and Result 

Let x be the center of the coated sphere (See Fig. ([I])). Inside the 
sphere, we ask that 

axApii = < |x - x| < r c 
a 2 Au = r c < |x — x| < r e , 

where a 1 and <r 2 are positive constants, together with the usual com- 
patibility conditions at the interfaces: 

u continuous across |x — x| = r c , (2.2) 

u = Ex\ at |x — x| = r e , (2-3) 
and the transmission conditions 

ffiii ■ |Vm| p ~ 2 Vm = o"2n ■ Vw, across |x — x| = r c , (2-4) 

and 

<T 2 n • Vw = (j^n • Vw, across |x — x| = r e . (2.5) 
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Using polar coordinates, we look for a solution u of ( 12. II) satisfying 
( i2.2p - O2.5p . with radial symmetry, of the form 

= a^r cos 9 for < r < r c , 

= £f cos 6* + a 2 r cos 9 for r c < r < r e , (2.6) 
= Er cos # for r = r e . 

where r = |x — x| and 9 measures the angle between the unit vector v 

in the direction of the applied field and x. 

We_ introduce the notation for the unit radial vector e r = n = 
x — x 

-, and denote by eg the unit vector perpendicular to e r . We 



|x — x| 

have v = cos#e r — sin^e^. 

Since (12.61) satisfies (12.11) . it is left to show that it satisfies the re- 
quired compatibility conditions at the interfaces, i.e., we require that 
the solution u of satisfies ([2T2]) - ([2~5j) . 

In what follows, we explain how the unknowns a±, a 2 , and b 2 and the 
effective conductivity a* are determined from (I2.2p . (12.31) . (12 .4p . and 

First, we look at the conditions u must satisfy when r = r c . From 
(J22D, we have 



o-i = a 2 + (2.7) 



r 



and from (12. 4p . we obtain 



cxi |ai| p 2 ai = a 2 ( a 2 - ^ ) . (2i 



We now look at the conditions that u must satisfy when r = r e . 
From (I2.3p . we have 

E = a 2 + %. (2.9) 



and from f ]2.5j) . we obtain 



^2 ( <*2 - ^ ) =<r>E. (2.10) 



Eliminating a 2 from ( 12 .7p and ( 12. 9p . we obtain 



«i = ^ + ^ (2 I1; 



Let us denote 
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and 



B 



2r 3 

~t + 1 



01 



+ 1. 



Observe that both A > and i? > are independent of r c and r e , 
i.e. they are defined only in terms of d\. 

Using (12. 8p and (12. 9p in (12.111) . we obtain the following identity 



0"! 



p-2 



°2 



which can be rewritten in terms of A and B as 



0"! 



E + A 



p-2 r 



E + A 



a 2 [E-B 



At this point, we consider the function 

f(x) =<ri\E + Ax\ p ~ 2 (E + Ax) - a 2 [E - Bx). 



(2.12) 



(2.13) 



Note that we obtain b 2 if we can prove that f(x) = has a (unique) 
solution. If that is the case, from (12. 9p and (12.111) . we can obtain a± 
and a 2 and from ( 12.101) . we can get an expression for a*. 
Let us study f(x). If E + Ax > 0, we have 

f(x) = a x {E + Axf- 1 - o 2 (E - Bx). 

Taking the derivative of the f(x), we have 

f\x) = Aatip -1)(E + Ax) p ~ 2 + a 2 B 

which is positive for all x, so the function f(x) is increasing. 
If E + Ax < 0, we have 



f(x) = —oi(—E - Axf- 1 - a 2 (E - Bx) 



and here 



f(x) = Aa^p - !)(-£ - Axf~ 2 + a 2 B 



is positive for all x so the function f(x) is also increasing in this case. 

Observe that as x — > — oo, the function f(x) approaches — oo and as 
x — > oo, the function f(x) approaches oo. 

Therefore, we conclude that f(x) has a unique solution xq. Moreover, 
observe that the coefficients of f(x) depend only on ax, E, <r 2 , 9%, and 
p, then 



x = — = K(a 1 ,a 2 ,E,9 1 ,p). 



(2.14) 
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Consequently, from ( 12. 14(1 we have b 2 = ^o r e which allows us to 
obtain a 2 and a\ from (12.91) and (12.111) . respectively. To obtain an 
expression for cr* we use (12 .9|) in (12. 101) as follows 

and from (12. 14j) . we get 

a* = ^(E-3K(a 1 ,a 2 ,E,9 1 ,p)) (2.15) 

Here, we would like to emphasize that (12.151) shows that cr* does 
not depend on r c or r e . Therefore, we realize that <r* is independent 
of scale. The rest of the coated spheres inserted into the matrix are 
chosen identical up to a scale factor, so that when inserted, the field is 
not disturbed. In this way, 9\ is the proportion of nonlinear material 
in the resulting assemblage of coated spheres and a* is its effective 
conductivity. 




Figure 3. Cross section of the assemblage of coated spheres. 



Remark 2.1. If p = 2, 

f{x) = a x (E + Ax) - o 2 [E - Bx) 
((T2 — o~i)E 

has a unique root Xq = — and in this case 

Aui + Ba 2 

0"2 /„ 3(<7 2 - (Ti)E 

0* = — \ E 



E \ Ao x + Ba 2 

W i( j 2 (ai -o- 2 ) . , 

3<j 2 + 6^2(^1 - 0-2) 
which is the Hashin-Shtrikman formula. 

3. Assemblage of Coated Disks: Statement of the 
Problem and Result 



Following the same method, we obtained the same results in the two- 
dimensional case of assemblages of coated disks. Consider a particular 
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coated disk of center x with core radius r c and exterior radius r e , 1 < 
r c < r e , subject to linear boundary conditions E • x = Exi, (where 
E = Ee 1 , with e 1 = (1,0) and x = (xi,a; 2 )) as a boundary condition 
to the exterior boundary of the disk, to find a solution u that solves 
(11.11) together with the usual compatibility conditions at the interfaces. 
The ratio of the inner disk radius to the outer disk radius is fixed, 
therefore the area fractions of both materials 9\ and # 2 are fixed: 

^2 



r?7r 



01 = -5- = -5, and 9 2 



rlix 



0i- 



(3.1) 



Accordingly, using polar coordinates, we look for a solution u of (12. ip 
satisfying (I2.2p - (l2.5p . with radial symmetry, of the form 



a\r cos 9 



for < r < r r 



— cos 9 + a 2 r cos 9 for 



r r < r < r P 



(3.2) 



with r, 9, e r and described in a similar way as in the previous section. 
The corresponding equations to (12. 7\i - (12. 101) are 



(ii = a 2 + 



<7i Oi 



7. IP-2 ~ 



ai = cr 2 a 2 



-E = a 2 + 



,2 ' 



02 a 2 - 



(3.3) 

(3.4) 

(3.5) 
(3.6) 



Eliminating a 2 from (13. 3 P and (13. 5p . we obtain 

b 2 



Let us denote 

> 

r 



C 



2 

I" 1 

c 



ai - 
1 



E + 



1 and D 



+ 1 



1 + 1. 

0i 



(3.7) 



Observe that C > 0, D > and both are independent of r c and r e 
(they are defined only in terms of #1). 
Using (13.4p and ( 13.51) in (13.71) . we obtain the following identity 

p-2 



0"! 



CT 2 




(3.8) 
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At this point, we consider the function 



g(x) = o-i 



E + Cx 



P * (E + Cx) - a 2 (E - Dx). 



We conclude that g(x) has a unique solution x Q . Moreover, observe 
that the coefficients of g(x) depend only on a x , a 2 , E, p, and 9%, then 

x = ^ = K(a 1 ,a 2 ,E,p,9 1 ). (3.9) 

' e 

Consequently, from (I3.9P we have b 2 = x§r 2 e which allows us to obtain 
a 2 and % from (13.51) and (13.71) . respectively, and 

a* = j(E-2K{a 1 ,a 2 ,E,p,9 1 f) (3.10) 

Remark 3.1. If p — 2, 

g{x) = a x (E + Cx) - a 2 (E - Dx) 

has a unique root x = -77" ~~~ — an d in this case 

C(Ti+L)(T2 

cr 2 / 2(a 2 - (7i)S 



E V C(7i + L>0 2 

2fli<r 2 (<ri -<7 2 ) 
= CT2+ 2a 2 + ^(a 1 -a 2 ) ^ 

which is the Hashin-Shtrikman formula. 

Remark 3.2. When p = 2, the effective conductivity a* of any isotropic 
composite of phases 1 and 2 (both linear in this case) satisfies the 
Hashin-Shtrikman bounds f |HS62al lMil02] ) 

d9 2 a 1 (a 2 - ax) d9 x a 2 {<J\ - a 2 ) 

01 + -; 7—, r > 0* > 02 + 



da 1 + 9 1 (a 2 - 01) da 2 + #2(01 - a 2 ) ' 

where d = 2, 3 is the dimensionality of the composite and it is assumed 
that the phases have been labeled so that a\ > a 2 . Thus, for p = 2, the 
coated sphere (12.161) and coated disk assemblages (13. lip with phase 
1 as core and phase 2 as coating attain the lower bound and with 
the phases inverted, attain the upper bound. They are examples of 
isotropic materials that, for fixed volume fractions 9\ and 9 2 — 1 — 
attain the minimum or maximum possible effective conductivity. 
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4. Discussion on cr* and x 

In this section, we discuss and analyze a* and xq and also their 
behavior with respect to p. 

First, observe that if we evaluate the function / (12 . 13[) at x = — -|, 
we obtain 



and if we evaluate / (12 . 1 3[) at x = -|, we obtain 



E + A 



P-2 

' £ + A 



=<^ + (ir)) Pl > a 

Therefore, we can conclude that a; (12.141) satisfies 

which implies i£ + Aco > and — > 0; and the effective con- 
ductivity a* (I2.15P satisfies 

0-2I <^<aJl + j 



or equivalently 



3(T2#1 3(72^1 

cr 2 < (J* < (Jo H 

O — (72 (72 



where A = j- - 1 > and B = |- + 1 > 0. If 6»i = 1 (the sphere 
is made only of nonlinear material), we have A = 0, B = 3, and 
/(x) = tri |£| p_2 (£) - a 2 (£ - 3s); from where x = (^-^\ P ~ 2 ) E an d 

cr* = cti |i?| p ~ 2 . If 6*i = 0, from above we obtain a* = a 2 . 

Let us study the case when E = 1. Observe that /(0) = ai — cr 2 > 0, 
therefore we can improve (14. ip and we obtain 

- < x < 0, (4.2) 

which implies < 1 + Ax < 1. Note that, in this case, 

dxo_ -a 1 (l + Ar )P- 1 ln(l + Ar ) >Q 
dp cti (p- 1)A(1 + Ax ) p - 2 + cr 2 B ~ ' 
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and it is equal to if A = which corresponds to the case when 9\ = 1. 
So Xo is an strictly increasing function of p, except when d\ — 1 (all 
nonlinear material), in which case it is a constant function xq = a2 ? ~° i 1 ■ 
If E — 1 and, for example, oi = 10, 02 = 1, and 9\ = 0.0, 0.2, 
0.4, 0.5, 0.6, 0.8, 1, the behavior of xq for different values of p can be 
observed in Fig H] and the values of xo for p = 1.1, 1.3, 1.6,2,2.7,4, 10 
can be found in Table 1. 




Figure 4. Plot of xq with respect to p for E = 1, o\ 
10, a 2 = 1, 0i = 0.00001, 0.2, 0.4, 0.5, 0.6, 0.8, 1 



0i/p 



1.1 



1.3 1.6 



2.7 



4 



10 



0.0 
0.2 
0.4 
0.5 
0.6 
0.8 
1 



-0.00 
-0.25 
-0.67 
-0.99 
-1.43 
-2.33 
-3 



-0.00 
-0.24 
-0.61 
-0.87 
-1.19 
-2.03 
-3 



-0.00 
-0.21 
-0.52 
-0.72 
-0.98 
-1.74 
-3 



-0.00 
-0.18 
-0.43 
-0.60 
-0.82 
-1.50 
-3 



-0.00 
-0.14 
-0.34 
-0.47 
-0.65 
-1.24 
-3 



-0.00 
-0.10 
-0.25 
-0.35 
-0.48 
-0.96 
-3 



-0.00 
-0.05 
-0.12 
-0.17 
-0.24 
-0.52 
-3 



Table 



Values of Xq, when E — 1, ai — 10, and 02 = 1. 



In this case, observe that when Q\ 
Since cr* = (x 2 (l — 3x ), we have 



1, x 



<9cx* 



1-10 



-3(7 2 -^ < 0, 



-3. 



dp dp 

and it is equal to if A = which corresponds to the case when 
#i = l. So cr* is an strictly decreasing function of p, except when 
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#i = l (all nonlinear material), in which case it is a constant function 
c* = <J\ (=10 with the values given previously). The behavior of a* for 
different values of p can be observed in Fig [5] and the values of a* for 
p = 1.1, 1.3, 1.6, 2, 2.7, 4, 10 can be found in Table 2. 

10-i 



9- 




2 ' 3 ' _ 4 ' 5 ' 6 7 ' 8 ' 9 "~~?0 
P 



Figure 5. Plot of a* with respect to p for E — 1, a\ — 
10, C7 2 = l,#i = 0.00001, 0.2, 0.4, 0.5, 0.6, 0.8, 1 
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1.72 


1.63 


1.53 


1.42 


1.30 


1.14 


0.4 


2.99 


2.84 


2.55 


2.29 


2.01 


1.74 


1.36 


0.5 


3.98 


3.62 


3.17 


2.80 


2.42 


2.05 


1.51 


0.6 


5.29 


4.57 


3.94 


3.45 


2.94 


2.45 


1.72 


0.8 


7.99 


7.09 


6.22 


5.50 


4.72 


3.89 


2.57 


1 


10 


10 


10 


10 


10 


10 


10 



Table 2. Values of a*, when E — 1, u\ — 10, and o 2 = 1- 

Also, as a function of 6i, the effective conductivity cr* satisfies 

fc* = ( ai(p-l)(l + Ax y- 2 + 2a 2 \ 

d6 1 " ° 2 Q\ Ui(p- 1)^(1 + Ax )p- 2 + Ba 2 ) 

Therefore cr* is an increasing function with respect to Q\. The plot 
of cr* varying with respect to the volume fraction 9\ in the case when 
E = 1, o\ = 10, and o-i — 1 can be observed in Fig [61 The blue plot 
corresponds to the case when p = 2 (Hashin-Shtrikman lower bound 
for different values of 9{). All the curves have fixed values a* = a 2 for 
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#i=0 and er* = a% for 9\ = 1. In the case when o\ = 10, and o 2 
then <r* = 1 for 9\ — and a* = 10 for 9 1 — 1 (see Fig [6]). 




Figure 6. Plot of a* with respect to 9\ for £7 = 1, ai — 
10, o- 2 = 1, p = 1.1,1.3,1.6,2,2.7,4,10. The blue plot 
corresponds to p = 2 (Hashin-Shtrikman curve), above 
this curve are the plots for p < 2 and below for p > 2 in 
decreasing and increasing order respectively. 

We now investigate the cases when < E < 1 and E > 1. For these 
cases, we have that 

dx -ai(E + AxoY" 1 \n(E + Ax ) 



Also, since a* 



dp ai{p - l)A(E + Ax Q )P~ 2 + a 2 B' 
- y-(E — 3x ), we have 



<9cr* _ g cr 2 dx 
dp E dp 

With respect to 6>i, we obtain 



dxo 



•TO 

01 



a x {p - 1)(E + Ax ) p - 2 + 2a 2 



(4.3) 



(4.4) 



(4.5) 



_<7i(p - 1)A(£ + Ai )P" 2 + a 2 5_ 

and the change of the effective conductivity cr* with respect to the 
volume fraction 9\ is given by 



do* —3o 2 x 



S0i 



£0? 



oi(p- !)(£ + Ax )P- 2 + 2cr 2 
_<7 X (p - 1)A(E + Ai )p- 2 + a 2 B_ ' 



(4.6) 



We start with < E < 1. In this case, since < + Aco < 1, 

we have from (14. 3ft that — — > 0, and we conclude xo is an strictly 

dp 
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2 


2.7 


4 


10 


0.0 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.2 


-0.17 


-0.17 


-0.15 


-0.12 


-0.09 


-0.05 


0.01 


0.4 


-0.47 


-0.45 


-0.38 


-0.30 


-0.21 


-0.12 


0.02 


0.5 


-0.70 


-0.65 


-0.54 


-0.42 


-0.29 


-0.17 


0.03 


0.6 


-1.04 


-0.92 


-0.74 


-0.57 


-0.40 


-0.22 


0.04 


0.8 


-2.17 


-1.73 


-1.35 


-1.05 


-0.73 


-0.41 


0.06 


1 


-2.99 


-2.75 


-2.45 


-2.10 


-1.58 


-0.91 


0.10 



Table 3. Values of xq, when E = 0.7, <j\ = 10, and a 2 = 1. 



increasing function of p. If E = 0.7, o\ = 10, a 2 = 1, and 6\ = 
0.0, 0.2, 0.4, 0.5, 0.6, 0.8, 1, the behavior of xq for different values of p can 
be observed in Fig[7[ and the values of xo for p = 1.1, 1.3, 1.6, 2, 2.7, 4, 10 
can be found in Table 3. 



Figure 7. Plot of x with respect to p for E = 0.7, 
a x = 10, a 2 = 1,6! = 0.00001, 0.2, 0.4, 0.5, 0.6, 0.8, 1 

By (14. 4p . we can also conclude that a* is an strictly decreasing func- 
tion of p. The behavior of a* for different values of p, in the case when 
E = 0.7, o"i = 10, and a 2 = 1, can be observed in Fig [8] and the values 
of cr* for p = 1.1, 1.3, 1.6, 2, 2.7, 4, 10 can be found in Table 4. 

The rate of change of the effective conductivity a* with respect to 
the volume fraction 6\ ( 14.6f) is determined by the sign of xq. If xq < 0, 
then a* is an strictly increasing function of 0\, and if xo > it is 
decreasing. We obtain have that <t* is a constant function of 9\ if 

— = — r . A plot of function of 6^ can be observed in 

a 2 Ep- 2 
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2 3 4 S 6 7 8 9 l'o 



P 

Figure 8. Plot of a* with respect to p for E = 0.7, 
Gl = 10, a 2 = 1,6! = 0.00001, 0.2, 0.4, 0.5, 0.6, 0.8, 1 
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2.26 


1.72 


0.89 


0.6 


5.47 


4.96 


4.15 


3.46 


2.70 


1.96 


0.85 


0.8 


10.3 


8.42 


6.79 


5.51 


4.13 


2.76 


0.74 


1 


13.8 


12.8 


11.5 


10 


7.78 


4.90 


0.58 



Table 4. Values of a*, when E = 0.7, o\ = 10, and a 2 = 1. 



Fig [9] for the values E = 0.7, <j\ = 10, and o 2 = 1. The blue plot 
corresponds to the case when p = 2 (Hashin-Shtrikman lower bound 
for different values of 6\). The curves have values cr* = a 2 = 1 for 
#i = and a* = a x \Ef 2 = 10(0.7) p " 2 for 6i = 1. Observe for 
example, for p = 1.1 and 9\ = 1, cr* = 10(0. 7)~ ' 9 ~ 13.8 and for 
p = 10, a* = 10(0.7) 8 ps 0.58. Observe in Fig M that for p = 10, 
the effective conductivity a* is decreasing with respect to 9%, this is 
obtained from the fact that xq > for p = 10, E = 0.7, a\ = 10, and 
er 2 = 1. 

When E > 1, we have — — > iff o\ > a 2 ( ) , otherwise 

dp V 1 - 0i / 

^ < 0. For example, if E = 2, a x = 10, a 2 = 1, and Q x = 0.0, 0.2, 
op 

0.4, 0.5, 0.6, 0.8, 1, the values of x for p = 1.1, 1.3, 1.6, 2, 2.7, 4, 
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Figure 9. Plot of a* with respect to 9 X for E = 0.7, 
oi = 10, a 2 = 1, p = 1.1,1.3,1.6,2,2.7,4,10. The blue 
plot corresponds to p = 2. 



10 can be found in Table 5. Observe that, for instance, for 9\ = 0.4, 



xq is increasing as p increases, but for 9\ 
increases. 



0.8, xq is decreasing as p 



Oi/P 


1.1 


1.3 


1.6 


2 


2.7 


4 


10 


0.0 


-0.00 


-0.00 


-0.00 


-0.00 


-0.00 


-0.00 


-0.00 


0.2 


-0.49 


-0.43 


-0.39 


-0.35 


-0.32 


-0.30 


-0.27 


0.4 


-1.14 


-1.01 


-0.92 


-0.86 


-0.80 


-0.76 


-0.70 


0.5 


-1.48 


-1.35 


-1.26 


-1.20 


-1.14 


-1.09 


-1.04 


0.6 


-1.79 


-1.73 


-1.67 


-1.64 


-1.59 


-1.56 


-1.52 


0.8 


-2.38 


-2.56 


-2.78 


-3. 


-3.25 


-3.5 


-3.81 


1 


-2.90 


-3.43 


-4.40 


-6. 


-10.1 


-26. 


-1710 



Table 5. Values of xq, when E = 2, o~i = 10, and <t 2 = 1. 



The values of a* for p — 1.1, 1.3, 1.6, 2, 2.7, 4, 10 can be found in Ta- 

/ 2^ 2 Q 

ble 6. By ( 14.41) . we have cr* it is increasing iff o\ < a 2 



otherwise 



dp 



1 



> 0. For example, for 9\ = 0.4, the value of a* decreases 



as p increases but for 9\ = 0.8, <r* increases as p increases. 

Observe that /(0) = axE^ 1 - a 2 E > if E > 1, therefore x < 
and, by (14. 6p . the effective conductivity cr* is an strictly increasing 
function of 9\. The effective conductivity function of the volume 

fraction 9\ for the case when E — 2, <j\ — 10, and a 2 — 1 can be 
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Bi/p 


1.1 


1.3 


1.6 


2 


2.7 


4 


10 


0.0 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


0.2 


1.73 


1.65 


1.58 


1.53 


1.48 


1.44 


1.40 


0.4 


2.71 


2.52 


2.38 


2.28 


2.20 


2.13 


2.05 


0.5 


3.22 


3.02 


2.89 


2.80 


2.71 


2.64 


2.56 


0.6 


3.68 


3.60 


3.50 


3.46 


3.38 


3.34 


3.28 


0.8 


4.57 


4.84 


5.15 


5.50 


5.90 


6.25 


6.71 


1 


5.35 


6.15 


7.60 


10 


16.2 


40.0 


2560 



Table 6. Values of <r*, when E — 2, a\ — 10, and oi = 1. 



observed in Fig [TUJ The blue plot corresponds to the case when p 
(Hashin-Shtrikman lower bound for different values of 9i). 




Figure 10. Plot of a* with respect to Q\ for E = 2, 
a x = 10, a 2 = 1, p = 1.1,1.3,1.6,2,2.7,4,10. The blue 
plot corresponds to p = 2. 



5. Final Remarks 

We considered the problem of constructing neutral inclusions from 
nonlinear materials. In particular, we studied the case of a coated 
sphere in which the core was nonlinear and the coating was a linear 
material. We showed that the coated sphere is equivalent to a sphere 
composed only of linear material of conductivity cr* f )2.15p . Thus the 
coated sphere is neutral in an environment formed by a linear mate- 
rial of conductivity a*. One could continue to add coated spheres of 
various sizes but with fixed volume fraction 1( without disturbing the 
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prescribed uniform applied electric field surrounding the inclusions, fill 
the entire space with an assemblage of these coated spheres and this 
configuration of nonlinear materials dissipates energy the same as a 
linear material with thermal conductivity cr*. We then studied the 
two-dimensional case of assemblages of coated disks. 

Ongoing work includes the generalization of this work to assemblages 
of coated ellipsoids, and to investigate if the coated spheres (disks) as- 
semblage constructed in this paper could be an optimal microstructure, 
in the sense that these microstructures attain bounds for the effective 
conductivity of a composited made of a linear matrix filled with a p- 
harmonic material. 
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